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H ' Abstract 

We investigate the planar solution of matrix models derived from various Cliern-Simons-matter theo- 
ries compatible with the planar limit. The saddle-point equations for most of such theories can be solved 
in a systematic way. A relation to Fuchsian systems play an important role in obtaining the planar 
resolvents. For those theories, the eigenvalue distribution is found to be confined in a bounded region 
even when the 't Hooft couplings become large. As a result, the vevs of Wilson loops are bounded in 
the large 't Hooft coupling limit. This implies that many of Chern-Simons-matter theories have quite 
different properties from ABJM theory. If the gauge group is of the form \J{Ni)ki x U(A^2)fe2i then the 
resolvents can be obtained in a more explicit form than in the general cases. 
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1 Introduction 



Supersymnietric Chern-Simons-niattcr theories have been intensively studied for several years. A recent 
motivation for studying these theories is to understand the dynamics of M2-branes. The worldvolume 
theory on M2-branes in the flat space-time was constructed in [I] [2] [3] [1] [5] whose applicability is, however, 
limited to the system with two M2-branes. This limitation on the number of M2-branes was lifted in 
[6] in which the background space-time is replaced with an orbifold. The worldvolume theory, known as 
ABJM theory, allows us to take the large N limit, where N is the number of M2-brancs. In this limit, one 
finds a relation to M-theory on AdS4 x S'^/Z^ [HI in a quite similar way to the well-studied AdS5/CFT4 
correspondence [7]. 

Later developments of the study on ABJM theory and other supersymmetric Chern-Simons-matter 
theories were enabled by the derivation of the localization formula for the partition functions and vcvs of 
some BPS operators of these theories [5]. For ABJM theory, the localization formula 
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has been studied in various ways. One of the most apparent way to tackle this formula would be to 
regard this integral as a partition function of a matrix model, and to apply the techniques developed 
in the context of matrix models. Indeed, such researches appeared soon after the derivation of (|l.ip 
[9] [To] . In [To] , a relation to a topological string theory was also found which turned out to provide a very 
effective method to derive various results [11] . The power of the relation to the topological string theory 
was also shown in fl2| which succeeded in deriving a surprisingly simple formula for the free energy of 
ABJM theory which includes all 1/A^ corrections. The same formula was also derived in [13] in a different 
way. In |13| . the partition function (jl.ip was related to the thermodynamical partition function of a 
non-interacting Fermi gas. Such a relation to a Fermi gas is familiar in the context of matrix models, but 
for the case of ABJM theory, the one-particle Hamiltonian is quite non-trivial. The relation to a Fermi 
gas can be extended |13j to various Chern-Simons-matter theories of quiver types originally studied in 
|14j . In fact, the localized partition function of a generic Chern-Simons-matter theory can be written 
as the partition function of a Fermi gas |15| . but it is interacting in general, and the analytic power of 
the reformulation seems to be limited. Another type of analysis was developed in |14j which restricts 
oneself to the leading behavior in a particular large N limit. (See also [9] for a related analysis.) It 
was shown [14] [16] [17] that this technique reproduces the volume of the internal manifold in the dual 
M-theory background. 

In [18], we investigated the large N solution of a particular family of theories, called Chern-Simons- 
adjoint theories. This family consists of A/" = 3 lJ{N)k Chern-Simons theories coupled to arbitrary 
numbers of adjoint hypermultiplets. For most of these theories, the Fermi gas approach [13] results in 
interacting theories. It was shown in [18] that if we restrict ourselves to the planar limit, the solution 
can be obtained explicitly in terms of theta functions. The planar resolvent is given as an integral of 
a combination of theta functions. Although the resulting expression does not look so simple, it was 
shown [18] that some properties of the eigenvalue distribution and the behavior of the vev (W) of a BPS 



Wilson loop as a function of the 't Hooft coupling t can be investigated in some detail. For example, 
(W) turned out to be finite even when t diverges. This is a quite different behavior from the one which 
would be expected if there would be a classical AdS dual and (W) would be given in terms of a minimal 
surface in the AdS space [H] [20] [21] . It would be natural to wonder whether this behavior is special to 
Chern-Simons-adjoint theories. To clarify the situation, it is desired to extend the analysis in [TH] to 
more general theories. 

In this paper, we discuss all A/" = 3 Chern-Simons-matter theories which are compatible with the 
planar limit. Since the constraint due to A/" = 3 supersymmetry is not so severe, there are huge number 
of such theories. The ability to take the planar limit requires that the representations of the matter 
multiplets should be small enough. It turns out that the analysis of a generic theory is reduced to the 
one of a related theory including only fundamental and bi-fundamental multiplets, under an assumption 
on eigenvalue distributions based on the numerical observations in |14] . The planar resolvent of the 
simpler theory can be written in terms of the solution of a system of differential equations unless the 
parameters of the theory satisfy an algebraic equation. Although the solution cannot be written in terms 
of simple functions, the properties of the solution can be investigated in enough detail. In order to relate 
the problem of solving saddle-point equations to a system of differential equations, the knowledge of the 
Riemann-Hilbert problem in the most basic setup is used. We mainly focus on the behavior of (W) for 
large t. If a Chern-Simons-matter theory has a simple gravity dual in some limit, then (W) should exhibit 
an exponential behavior as t becomes large. It would be more interesting to investigate, in addition, the 
behavior of the free energy. This will be hopefully discussed elsewhere. 

It turns out that any Chern-Simons-matter theory, which is solvable by the technique developed 
in this paper, exhibits the behavior similar to Chern-Simons-adjoint theories [18]. This suggests that 
Chern-Simons-matter theories with simple AdS gravity duals would be quite rare even though there are 
huge number of superconformal Chern-Simons-matter theories. A condition is found which shows when a 
Chern-Simons-matter theory might have more interesting properties. Such "exceptional" theories include 
ABJM theory and all other theories with known gravity duals. It would be interesting to compare this 
situation with the one in AdS5/CFT4 correspondence. In four-dimensions, superconformal field theories 
are already rare, but many of them have gravity duals. 

This paper is organized as follows. In section [21 we define a family of theories which is investigated 
in this paper. Starting with the localization formula for the partition function, we derive the saddle- 
point equations. Assuming some symmetry for the eigenvalue distributions, the problem of solving 
the saddle-point equations for a generic theory is reduced to the one for a related theory including 
only fundamental and bi-fundamental multiplets. Unless the theory is "exceptional" , the saddle-point 
equations can be further simplified. In section [S] we focus on theories for which the gauge group is of 
the form U(iVi)fej x U(A^2)fc2- The saddle-point equations for these theories can be solved explicitly in 
terms of theta functions, as in |18| . The condition for the existence of the solution is investigated, and 
this turns out to be the condition for a divergenct behavior of the 't Hooft couplings. In section JH we 
solve the saddle-point equations for more general theories. It turns out that the technique used in section 
[3] cannot be straightforwardly extended to the general cases. Instead, the saddle-point equations can be 
solved by relating them to the solution of a Fuchsian system, a set of first-order differential equations. 
The relation to the Fuchsian system is given by the known results on the Riemann-Hilbert problem. 



Using the structure of the solution of the Fuchsian system, it is possible to extract some properties of 
the eigenvalue distributions. A relation between the divergent behavior of the 't Hooft couplings and the 
isomonodroniic deformation of the Fuchsian system is pointed out. Section [5] is devoted to discussion. 
Appendices contain some technical details. 



2 Planar limit of J\f = 3 Chern-Simons-matter theories 



A family C of theories discussed in this paper consists of A/" = 3 Chern-Simons-matter theories which are 
compatible with the planar limit. In this section, we describe the family C, and for each theory in C, 
derive the equations which determine the saddle-point of the localized partition function. Some physical 
quantities can be obtained from the solution of the saddle-point equations. It will turn out that the study 
of a generic theory in C in the planar limit is essentially reduced to the study of a theory in a sub-family 
Co in which the theories consist of only fundamental and bi-fundamental multiplets. 



2.1 Saddle-point equations for theories in C 

A theory in C is an A/" = 3 Chern-Simons theory coupled to various matter multiplets. The gauge group 
G of the theory is a direct product of several unitary groups, 

rig 

G = nU(iV„)fc„, (2.1) 

a=l 

where ka is the Chern-Simons level for the gauge group factor \J{Na)- The gauge fields with their 
superpartners form A/" = 4 vector multiplets, and the matter fields form A/" = 4 hypermultiplets. The La- 
grangians of Chern-Simons-matter theories with A/" > 4 supersymmetry were constructed in [22] [23] [6] [24] . 
For a generic choice of G and matter multiplets, a Chern-Simons-matter theory possesses at most Af ~ 3 
supersymmetry since the Chern-Simons terms break the Af = 4 supersymmetry partially. For an A/" = 3 
theory, the superpotential is uniquely determined by the symmetry. Therefore, the Lagrangian of an 
A/" = 3 theory is uniquely specified by the choice of the gauge group G including the Chern-Simons levels 
{ka} and the representations of G to which the matter multiplets belong. One can find more details on 
supersymmetric Chern-Simons-matter theories in |25j . 

Note that, for some choices of the data mentioned above, it is argued that the corresponding A/" = 3 
Chern-Simons-matter theories do not exist as consistent quantum field theories. An example is ABJ 
theory [26], regarded as an A/" = 3 theory, with the gauge group U(A^i)fe x U(A^2)-fc where the parameters 
satisfy \Ni ~ N2\ > k. We assume that such theories are excluded from C. It should be noted that 
the criterion for selecting the data which provide us with a consistent theory does not seem to be fully 
investigated. This issue will be discussed elsewhere. 

As mentioned above, each theory in C is further restricted such that it is compatible with the planar 
limit. The possible representations of G allowed for theories in C are the followings. 
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Figure 1: An example of a diagram with two nodes. Various edges correspond to matters in various 
representations. 

• a fundamental multiplet (f) in Na ® Na of \J{Na), 

• an adjoint multiplet (ad) in adj adj of \J{Na), 

• a symmetric multiplet (s) in sym sym of \J{Na), 

• an anti-symmetric multiplet (as) in asym ® asym of \J{Na), 

• a double-fundamental multiplet (ff) in {Na «) Nb) ® (iVa® iVb) of U(iVa) x U(A^6), 

• a bi-fundamental multiplet (bf) in {Na <»1^) © {7^(BNb) of U(iVa) x U(iVb). 

Note that an A/" = 4 hypermultiplet consists of two M = 2 chiral multiplets in the representations R and 
i?. 

Each theory in C is associated with a diagram F which consists of n^ nodes corresponding to U(Afa) 
factors and a number of colored edges each of which corresponds to one matter multiplet. An edge with 
a free end corresponds to f. An edge connected to the same node corresponds to either ad, s or as which 
are distinguished by a suitable choice of colors. An edge connected to two different nodes corresponds 
to either ff or bf. An example of such a diagram is shown in Figure [TJ Since the hypermultiplets are 
non-chiral, there is no natural rule to assign directions to the edges. In other words, F is not a quiver 
diagram. In the following, a theory in C associated to a diagram F will be often called a theory F. 

For a given theory F, one may calculate the sphere partition function Z{T)^ that is, the partition 
function of the theory F defined on the Euclidean three-sphere S^ . Using the supersymmetric localization 
[8], ZiV) can be written in terms of a finite-dimensional Riemann integral. The integral is taken over the 
eigenvalues {u^ } (i^ = 1, • • • , A^a) of Ug matrices {cr""} where each cr" is an Na x Na Hermitian matrix. 
The explicit formula for ^(F), with the overall numerical constant omitted, is as follows. 



where S'j^.pciu] comes from the classical action, S^[u\ is the contributions from the a-th vector multiplet 
and 
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is a sum of contributions from the matter multiplets. Here n^ and n'^ are the numbers of the matter 
multiplets in the representation R coupled to \]{Na) and \]{Na) x U(7V(,), respectively. The details of 
those terms are as follows, 
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In the following, we replace S'°j(,(,[m] with 
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and assume that all ka are real and positive. This is sufficient to make the integral (|2.2p well-defined 
for any choice of matter multiplets. The partition function of the original theory will be recovered via a 
suitable analytic continuation of {ka\ from real positive values to purely imaginary values. 

Consider a limit in which {Na\ and {ka\ are all proportional to a common number, say k, and k 
becomes infinitely large. The proportionality constants are denoted as 
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(2.5) 



where {i"} are the 't Hooft couplings of the Chern-Simons-matter theory. If there are a number of fs, 
{n^ } are also assumed to be proportional to k. As is familiar in the context of matrix models, the integral 
Z{T) is dominated by a particular configuration of the eigenvalues {u° } in this limit. The configuration 
is determined by the saddle-point equations 
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Let u denote a solution of these equations. The free energy F(r) of the theory F in the planar limit 
is given as 

"g "g 

F(r) = g^^rceN+E^vN+^mN- (2-7) 
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The leading term of F{T) in the planar limit is proportional to k^, and the coefficient is a function of the 
't Hooft couplings {t°-}. Some other observables consistent with the localization can be also determined 
in terms of u. For example, the expectation value of a supersymmetric Wilson loop for a \J{Na) factor is 
given as 

<^°) = ^^E'^^pK)- (2.8) 
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The theory F may have interesting properties, for example, the existence of a classical gravity dual. 
Typically, such a dual description appears when (some of) {t°} become large. Since the solution u is 
determined once {t°} are given (it may not be unique, see section |31), one can regard (VF°) as a function 
of {t"}. If {[¥"■) would be given in terms of a minimal surface in the dual geometry (possibly an AdS 
space), then {W"") would grow exponentially as {i"} grow. Therefore, the behavior of (VF") when {t°} 
diverge indicates whether the theory F is interesting in the above sense. 



2.2 Planar relations in C 

The family C defined so far consists of various kinds of Chern-Simons-matter theories. The variety looks 
so rich that one might think one needs to investigate those theories one by one. Fortunately, it turns 
out that they can be investigated in a rather systematic way under an assumption on the eigenvalue 
distributions. 

In [14j . the eigenvalue distributions were numerically determined for some matrix models of Chern- 
Simons-matter theories. The plots in [14) shows that the distributions are invariant under the reflection 

< -> -<• (2-9) 

In fact, this reflection, performed simultaneously on all eigenvalues, is a symmetry of the saddle-point 
equations (j2.6p . and therefore, it would be natural to assume that the solution should be also symmetric. 
In the following, we assume that 

{u'i,--- ,u%J = {„«-,. ..,_«^J (2.10) 

holds as sets for any a = 1, ■ ■ ■ ,ng. Under this assumption, it turns out that different Chern-Simons- 
matter theories share the same saddle-point equations, and therefore the same eigenvalue distributions. 
For example, the contribution to (|2.6p from ff, 



_^::^tanh^^^, (2.11) 
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coincides with the contribution from bf, 
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with n^f = n^. Therefore, as long as the planar solution is concerned, one may replace ff with the same 
number of bf and obtain the same eigenvalue distributions. In this sense, it is not necessary to discuss 
any theory F including ff. 
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Figure 2: The procedure for replacing ad with bf. 



By the same token, the matters s and as, whose contributions are 



2 



\^ tanh ■ 



■ tanh u? 






Oa 



can be replaced with ad, which contributes as 
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2 



(2.13) 
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(2.14) 



with nl^ = in^ 



up to "local terms" 



■ tanh 1 



(2.15) 

In fact, it can be shown that these "local terms" can be ignored in the planar limit. If n" and n°g are of 
order A;°, then this term is smaller than the other terms in (|2.6p by k~^ . On the other hand, if either n^ 
or n°g is of order fc, then the term (|2.14p becomes larger than the other terms by fc, and the equations 
p.6p do not have a non-trivial solution. 

There is a further simplification in which ad can be replaced with bf by changing the diagram F. This 
possibility can be anticipated from the similarity between ()2.12|) and ()2.14p . Suppose that, in a theory 
F, there is an a-th node with edges for ad. The replacement goes as follows. One eliminates the edges 
in F corresponding to ad and obtains a diagram F'. Then one duplicates F' and connects the a-th node 
in the original F' with the same node in its mirror F' by the edges corresponding to bf. The number of 
edges for bf is the twice the number of ad which originally existed. See Figure [2J Denote the resulting 
diagram as F" . The saddle-point equations (|2.6p for the theory F" are no longer equivalent to the original 
ones. However, they are equivalent to the original equations for the theory F if the mirror pairs of nodes 
in F" are assumed to have the same eigenvalue distributions. In this way, the planar solution of the 
original theory F can be obtained from the solution of the enlarged theory F" by suitably constraining 
the eigenvalue distributions. Repeating this procedure, one can find a diagram F'" including only edges 
for f and bf from which one obtains the planar solution of the original theory F including ad. 



Now, the problem of solving saddle-point equations (|2.6p for a generic theory in C has been reduced 
to solving 



h^u'i +-tanh^ = y coth^^^^-VV — tanh^^^^^^, (2.16) 
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where n°,ri,°'' are integers. Note that n'^,Ti°-^ are always non- negative. These are the saddle-point equa- 
tions for a theory F including only edges for f and bf, that is, the equations for a theory in Cq. Note that 
the analysis below shows that n" , n"''' can take any values other than integers as long as the equations 
(|2.16p are concerned. In fact, there exists a solution even when they are complex numbers. 



2.3 Equations in terms of resolvents 

To study the equations p.l6|) further, it is convenient to introduce new variables 

zl := exp«). (2.17) 

The reflection symmetry (j2.9p of uf is translated into the symmetry under the inversion 

< ^ «)-'• (2.18) 

In terms of zf , the equations (|2.16p can be written as 

^logz- +!Ll!.^^ = y ti^±l2^_y^y'_i^^_ (2.19) 

2^^ " 2<^ + l Z^^zi-zi Z^2J^zl+zl ^ ' 

A standard strategy to deal with this kind of equations is to rewrite them in terms of resolvents 
defined as 
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The 't Hooft coupling i" can be easily read off from v°'{z) as 

r = -w"(0) = u"(cx3). (2.21) 

The vev {W"") can be also obtained from the expansion of v°-{z) for large z, 

z;"(z)= t''[l + 2{W'')z-^ + 0{z-^)]. (2.22) 

Note that the symmetry under (|2.18p implies 

v^iz-^) = -v'^iz). (2.23) 

In the planar limit, it is assumed that, for each a, the eigenvalues {zf } condense to form a continuous 
distribution whose support is an interval [pa, Qa], and u°(z) becomes a holomorphic function on ¥^\[pa, Qa] 
having a branch cut on [pa , Qa] ■ Note that the symmetry under the inversion (j2.18p implies 

PaQa = 1 (2.24) 




Figure 3: Two limits y^ to approach ya- 



for each a. Wc further assume that f"(z) is finite at the branch points. This is based on the integral 
expression of v°'{z) 



v\z) 



dx p''{x) 



Z + X 



(2.25) 



where p'^{x) is assumed to be continuous in the planar limit. As long as p°(x) decays near x = Pa a.s a 
positive power oi x — Pa, then v°-(z) is finite at z = Pa- 

The right-hand side of (|2.19p can be written in terms of v°-{z) as 



«"(y.+) + «°(ya-)-E""'"'(-y") 



(2.26) 



where ya S [pa, Qa]- The quantities y^ are defined to be ya ± iQ if the interval [pa, qa\ lies on the real axis. 
Otherwise, y^ are defined as limits as depicted in Figure |3l The expression (|2.26p is, however, not useful 
to be dealt with since it is "non-local", that is, the functions in (|2.26p depend on both ya and —ya- A 
better expression can be obtained by introducing another set of resolvents 



1 



^ ' Na^ z + z'^ 

3a. J" 



Note that 



v''{z) = v''{-z) 



holds. Then, obviously (|2.26p becomes 



^°(2/a+)+«°(ya)-E""'^''(2^-) 



(2.27) 



(2.28) 



(2.29) 



which looks "local" . This is also equivalent to another "local" expression 

\ i-{-y+) + V-{-y-)^Y.n'^^v\-ya) 



(2.30) 



These expressions suggest that the equations (j2.19p can be written in "local" forms if one can choose 
either v"'{z) or v°-{z) for each node consistently such that any edge in the diagram F connects a node for 
which v°'{z) is chosen and another node for which v°'{z) is chosen. This choice is possible if and only if 
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Figure 4: The bipartite double cover of a triangle. 

the diagram F is bipartite. In this case, one may choose v'^{z) for white nodes and u''(z) for black nodes. 
Then the equations (|2.19p can be written as 



2K°l0gJ/a 

2K''\og{-yb) + 



n-ya 



k Va + 1 
k -yt + 1 



«"(2/^) + "°(ya")- 
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(2.31) 


*'(%+)+«'(%)- 


-)^n'V{y,), 


(2.32) 



where a runs over white nodes and b runs over black nodes, and y^ E [^Ib, ^Pb]- Since, for each a, either 
v°-{z) or v°-{z) appears in the above equations, one may forget the constraints (|2.28p . 

In the case when T is not bipartite, for example when F is a triangle, one may obtain a bipartite 
diagram from F by duplicating the nodes and draw edges such that each edge connects a node in the 
original diagram and another node in the mirror. See Figure 21 The resulting diagram F, known as the 
bipartite double cover of F, is always bipartite by construction. The saddle-point equations for F are 
nothing but (|2.3ip(|2.32p with the indices a, b both run over all nodes of F. One may obtain the resolvents 
of the original theory from the solution of the equations (|2.31l) (|2.32p by imposing the constraint (|2.28p . 

In the following, the diagram F is always assumed to be bipartite. The equations (|2.3ip(|2.32p can be 
written simply as 

2«°log(6(a)yJ + ^ ^|°j^°~| = ."(y+) + .«(ya-)-E"°'^'(2^'')' (^.33) 



b 



where the branch cut of w°(z) is on the interval [pa, Qa] and 



e(a) = I +'' ('^ " " "^^*° "°^^') (2.34) 

I —1, (a is a black node.) 



2.4 Integral ansatz 

The equations p.33p look like the ones encountered in the context of matrix models, except for the 
logarithmic functions in the left-hand side. They have a branch cut, which would make the procedure of 

11 



solving the equations, based on the analyticity, more difficult than the usual cases where only polynomial 
functions appear. 

It was found in TBI that the following ansatz 



v'^iz) = / d^z;°(z,0 (2.35) 

Jc 

for the resolvents is quite useful to deal with the logarithmic functions. Here the integration contour C 
starts from the origin of the complex plane C and reaches infinity, avoiding the branch cut of v°'{z). 

The equations (|2.33p can be reduced to a set of homogeneous equations using (|2.35p as follows. Let 
us first consider the case n° = 0. Suppose that a set of functions {u''(z,^)} satisfies the equations 

r{ya.i) = ^'''(y+,e) + «'^(2/a,e)-E""'«'(y-?) (2-36) 

b 

where 

ri^,0 = -Tl^-F^ (2.37) 

e{a)z-^ C-1 

are rational functions of z. Integrating both sides of the equations (j2.36p . one recovers the original 

equations (|2.33p . 

The integrand v°'{z,^) is assumed to be holomorphic on P^\[p",(7°], to have a branch cut on [p°-,q°-] 
for all values of <^ and to be finite at the branch points. As long as the integral p.35p is well-defined, the 
function v"-{z) obtained from the solution of (|2.36p is also holomorphic on P^\[p°,(7''] having a branch 
cut on [p°-, q"-], and finite at the branch points. 

To eliminate the rational function f°'{z,^), let us define a complex function a;"(z,^) such that 

v%z,S,) = r''(z,e)+c.'^(z,0 (2.38) 

holds for some rational function r'^(z,<^). If the following set of linear equations 

r(z,0 = ^(2<5«^-7i''V(^,0 (2.39) 

b 

has the solution, then the equations (|2.36p are reduced to the following homogeneous equations 

c.«(y+,0+^"(j/a",0-E""''^'(y«) = 0- (2-40) 

b 

The function a;°(z,^) must have appropriate poles which cancel those of r°(z,^) such that w"(z,^) does 
not have such poles. Note that the introduction of poles is much easier to be dealt with than that of 
log-cut, which would be necessary at the level of w°(z). This is an advantage of the integral ansatz (|2.35p . 
In the case n" > 0, the reduction may proceed in two steps. Since the terms in (|2.36p proportional to 
n° are rational, they can be eliminated by assuming 

v''{z) = rl{z) + v1{z) (2.41) 

where {?'i(^)} are rational functions satisfying the equations similar to (|2.39p . Then, employing the 
integral ansatz, the remaining logarithmic terms can be eliminated as above. 
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If the equations (|2.40p could be solved, then the 't Hooft coupling t°- could be obtained from u!°'{z, ^) 
as 

r = - /der(0,O+^"(0,O]- (2.42) 

Jc 

This equation determines t° as a function of the positions {pa, Qa} of the branch cuts. The analytic 

continuation of t° can be done by changing the configuration of the branch cuts. 

Note that the parameters which are given by specifying a diagram F are the 't Hooft couplings {i°}, 
while the positions {pa,qa} of the branch cuts would be determined according to the relations p.42|) . 
Later, it will turn out that the relations ()2.42p for all a do not always define a one-to-one map between 
the two parameter spaces, even when PaQa = 1 is imposed for any a. This implies that the saddle-point 
is not always unique for a given set of {i"}. In the case there are multiple saddle-points, it is necessary 
to compare the values of the free energy at those saddle-points to correctly extract physical information. 
This issue will not be investigated further in this paper. 

Now, the problem is reduced to solving the above homogeneous equations (|2.40p . In the next section, 
we investigate the simplest non-trivial case in which the diagram T consists of two nodes. Theories with 
more nodes will be discussed in section |4l 

It was found above that if 

det(2<5"'' - 71'''') = (2.43) 

is satisfied for a theory F, then the technique which will be developed in later sections cannot be applied 
to the theory. Indeed, since the left-hand side of p.39p is a set of functions, they cannot be always 
in the image of a degenerate constant matrix 2(5"'' — n"-''. We call a theory F in Cq satisfying (|2.43p 
a degenerate theory. Interestingly, many Chern-Simons-matter theories studied so far are degenerate. 
The examples include ABJM theory [6] [26] and its orbifolds [27], circular quivers [28], necklace quivers 
|14j and more general quivers |29j | 30j . These theories have been already solved by various techniques. 
Especially, some of them are investigated in detail using the Fermi gas approach [13] . See also [31] [32] [33] 
for recent achievements on instanton effects in this direction. It seems that those techniques are not so 
much useful to generic non-degenerate theories in Co, and therefore, the technique in this paper would be 
complementary to those known techniques. It is interesting to classify the solutions of p.43p and clarify 
what kinds of theories are included in the class of degenerate theories. The results obtained below for 
non-degenerate theories in Cq seem to suggest that the theories which have their simple gravity duals 
would be always degenerate. 

An example of an infinite family of solutions of (|2.43p which do not correspond to the theories men- 
tioned above is, for Ug = 4, 

^12 ^ „2i ^ ^34 ^ ^43 ^ 21, n23 ^ n^^ = 2^2-2 (2.44) 

and other components of n"'' are zero, for any positive integer /. It is interesting to clarify whether the 
theory corresponding to this data has interesting properties as other theories mentioned above. 
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3 Theories F with two nodes 

In this section, we discuss a simple class of theories in Cq whose diagrams contain only two nodes, that is, 
Ug = 2. Obviously, the corresponding diagrams are bipartite, and therefore, the saddle-point equations 
(|2.6|) can be written as (|2.33p . The equations (|2.33p can be reduced to the homogeneous ones (|2.40p if 
and only if n :~ n^^ does not satisfy 

4-71^ = 0. (3.1) 

The degenerate theories with Ug = 2 include, for example, ABJM theory [5], ABJ theory [55], ABJM 
theory with flavors [34] and GT theory [35]. The planar solution of ABJ(M) theory was obtained in [TO] . 
flavored ABJM theory was solved in [34] and GT theory was investigated in [36] . In the following, we will 
focus on the remaining cases n ^2. The technique for solving (|2.40p in this section is a generalization of 
the one developed in [57] } 38 ) . 

As explained in the previous section, the study of this class also gives us the planar solution of any 
theory T with a single node to which various edges corresponding to f, ad, s and as arc connected. 

The equations (|2.40p for Ug = 2 are 

u\y+)+uj\y^)~noj^iyi) = 0, (3.2) 

io\y+)+LU^{y^)-nio\y2) = 0, (3.3) 

where j/i G [pi,qi], y2 G [j'2,92]- For a while, the dependence on ^ will be ignored. The configuration of 
these intervals can be arbitrary in the complex plane C In fact, it is enough to discuss the case where 
the intervals arc symmetric, that is. 

Pi = -92, qi = -P2- (3.4) 

Indeed, there is an SL(2,C) transformation which maps the four points (pi, gi,p2,<?2) to {—q,—p,p,q). 
Explicitly, it is given as ipi o ip2 where 

I ■. 2; + ^/a , - q\- q2Z-vi 
ipi\z) := c ^, (p2{zj ■■= , (3.5) 



91 - 


- q2 z- 


-Pi 


91 - 


-pi z- 


-92 


y/a- 


1 





and 

a := (/32(P2), c^ := ^^7^7^- (3-6) 

The constant c is chosen such that 

pq = 1 (3.7) 

is satisfied. In the following, the intervals are assumed to be symmetric, satisfying the condition p.7p . 
The solution for a generic configuration can be obtained through the pull-back by (pi o (p2. 

By definition, uj^{z) and w^(z) do not have a branch cut on [p,q] and [—q,—p], respectively. This is 
expressed as 

u;\y+) = oj\y^), co\y+) = oj'{y^). (3.8) 

It is convenient to introduce a vector-valued function 

uj{z) := {uj\z),uj^{z)) . (3.9) 
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Figure 5: The red curve is the ^-cycle, and the blue eurve which passes through branch cuts is the 
i3-cycles. 



In terms of a;(z), the equations p.2pp.3pp.8p can be written as 

uj{y+) == uj{y-)Mi, uj{y+) = uj{y-)M2, 
where 



Ml 



Note that these matrices satisfy 



where / is the unit matrix. 



-1 
n 1 



, M2 



1 n 
-1 



Ml = Ml = I, 



(3.10) 



(3.11) 



(3.12) 



3.1 General solution 

As in |18| . it is convenient to change the variable. The new variable u has a geometric meaning |39| of 
an elliptic curve E defined by 

y^ = {x^-p'){x^-q% (3.13) 



Let A be the holomorphic 1-form on E normalized such that 



A = -r, / A = 1 

A JB 



(3.14) 



with Iin(r) > 0. The A-cycle and the S-cycle are depicted in Figure[51 This unconventional choice is due 
to the fact that we are interested in the large \t\ limit in the above definition where the 't Hooft couplings 
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Figure 6: The Abel-Jacobi map maps the whole z-plane with branch cuts into the parallelogram in the 
u-plane. The upper- and the lower-edges of the parallelogram are identified. The upper-half and the 
lower half of the left edge are the images of the same curve on C, lying on the different Riemann sheets. 
Similar for the right edge. 

are expected to diverge. In the usual convention, this limit correspond to r — > 0, and we have performed 
the iS-transformation from the beginning. In the opposite limit where the branch cuts shrink to points, 
one may perform the S'-transformation, or a perturbative calculation developed for ABJM theory in [9] 
can be adapted to provide a perturbative solution. 

The coordinate transformation to the new variable u is defined by the Abel-Jacobi map 



u(z) := / A 
Jp 

with the base point p. The image of the whole z-planc is a parallelogram 

K := <! a + ^r e C I < a < i, |/3| < i 



(3.15) 



(3.16) 



The horizontal edges of E^ are identified. Roughly speaking, the left edge corresponds to the union of 
the upper- and the lower-side of the branch cut [p,q\, and the right edge corresponds to [— g, —p]. See 
Figure |6l Some properties of the map are shown in Appendix \^ 



One can show that, if a vector- valued function J7(u) defined on the w-plane satisfies the following 
conditions 

n{u + i) = n{u)M2Mi, n{u + T) = n{u), n{-u) = vl{u)M2, (3.17) 



for any m e C, then uj{z) given by the pull-back, 

w(z) = r2(u(z)), 

is a solution of (|3.10l) . 



(3.18) 
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The equations p.l7|) can be solved easily by diagonalizing the matrix AI2M1. This can be done by a 
similarity transformation defined by a matrix 



S 

which is non-singular for 71 7^ 2. Here v is related to n as 

n = — 2cos7ri^. 
The transformed function ris(u) '■— ^(u)S satisfies, in terms of components, 

Using the identities of the theta functions, one can show that the following function 
G(m) := e^ 



2^^uu ^i {U-Uo + ^T, r)^i (u - Uq + i + |t, t) 



'di{u-uo,T)'di{u-uo + ^,t) 
where ug ^ C can be arbitrary, satisfies the equations (|3.2ip . Suppose that 5^5(w) is given as 

n],{u) = R{u)G{u). 



(3.19) 



(3.20) 



(3.21) 
(3.22) 



(3.23) 



(3.24) 



Then, the equations p.2ip imply that R(u) is an elliptic function on E. ri|(w) is given in terms of R(u) 
and G{u) by dSSS)- Exphcitly, 



nl{u) = -R{-u)G [-U + 



(3.25) 



Note that there is a simpler choice for G{u). A different G{u) just results in a different R{u). 



Now we recover the ^-dependence. It is obvious that the dependence of w(z, ^) on ^ is encoded in the 
elliptic function R{u,^) only. The functional form of R{u,£^) is constrained since w(z,^) must have the 
appropriate pole structure. 

For simplicity, we consider the case n" = first. The general case will be discussed in subsection 
3.51 In the case Ug = 2, the vector- valued function r(z,^) := (r^(z,^),r^(z,^)), whose components were 
introduced in (|2.38p . is given in terms of f{z,^) := (/^(z,^), /^(z,^)) as 



riz,0 



4-71- 



r/(^,0 



The function r(z,^) has poles at z = ±^ with the residues 

Res_^7'(z,^) 



Res+^r{z,0 = ^{n,2). 



A- n? 
4- n 



2 n 
n 2 



(2,n), 



(3.26) 

(3.27) 
(3.28) 
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Therefore, uj(z,^) ~ fl{u{z),^) must have poles with residues which have opposite signs to the above 
ones. In terms of ujs{z,^) := fls{u{z),^)^ the required residues are 

Res.^cusizA) = -^^^^(e-^^'^-e^-'''), (3.29) 

Res+^us{z,0 = -^^^^(e^"^-e-^-'''). (3.30) 

In the following, we abbreviate R{u{z),£,) as R{z,£_) etc. Since ujs{z,£,) must have poles at z = ±^, 
i?(z,^) can be written as 



^(^'^) = Z2^7^M^:0 (3.31) 



1 

with another elliptic function Ro{z,^) on E. To write down the most general form of Ro{z,^) explicitly, 
it is convenient to choose uq in G{z) as 

uq = Uoo := m(oo). (3.32) 

For this choice, the function 



F{z,0 ■■= -^ZTT^Giz) (3.33) 



1 

Z2~^2 

has simple poles aX z ~ ±^, simple zeros aX z — z^, oo, — Zjy, oo, where the tilded points lie on the second 
Riemann sheet. Here z^ is defined such that 

u{zi,) ^ Uoo- -zT. (3.34) 

The elliptic function Rq{z,^) may have poles which can be canceled by the zeros of F{z,^). In terms of 
the divisor, 

iRo{z)) > -(z,) - (oo) - (^,) - (55). (3.35) 

Recall that w(z,^) is assumed to be finite at branch points. 

The Riemann-Roch theorem implies that there are four such functions linearly independent on C. 
Therefore, Ro{z,^) can be written as 

4 

Ro{z,0 = E^*(^)/»(^) (3.36) 

for some basis functions fi{z). 

In summary, it has been shown that the general form of tJs{z, ^) is 

4(^:0 = ^^Giz)[c,i0hiz) + C2{0f2{z) + C3{0f3iz) + cmh{^)]. (3-37) 

^1(^,0 = -^^^G'(-z)[ci(e)/i(iO + c2(e)/2(i) + c3(e)/3(2)+c4(e)/4(i)]- (3.38) 

The four parameters Ci(^) can be determined by the four conditions for the residues, two for ujg{z, ^) and 
two for aj|(z,^). Therefore, the vector-valued resolvent v{z) obtained from u!s{z,£,) can be completely 
determined by the data given so far. 

It is important to note that the above analysis can be used to show that the (vector- valued) solution 
v{z) of the saddle-point equations (|2.33p is unique. Suppose that there are two solutions v(z) and v'{z). 
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The difference 5v{z) :— v{z) — v'{z) is a solution of the homogeneous equations (|2.40p . Therefore, the 
function dvg{z), which is the first component of dv{z)S, can be written as 

Svliz) = R{z)G{z). (3.39) 

In this case, R{z) does not have poles at z = ±^. It must have zeros sd z = 00,00 to cancel the poles 
of G{z). Since R{z) is an elliptic function, there must be at least two poles. However, the only allowed 
positions of the poles of R{z) are at z = Zj,, —Zi, since Svs{z) should not have any poles. The divisor of 
R{z) is therefore 

(i?(z)) = -(z,)-(^) + (oo) + (55). (3.40) 

Another constraint for elliptic functions is that the u-coordinates of these zeros and poles must satisfy 

- m(zi.) - u(^) + u(oo) + u(oo) e A, (3.41) 

where 

A := Z + Zr. (3.42) 

This then implies 

i/r e A. (3.43) 

Therefore, for generic values of v and r, the only allowed solution is 5v{z) = 0, implying the uniqueness of 
the solution of p.33p . It will be shown below that if p.43p is satisfied, then there is actually no solution. 



3.2 The solution in detail 

It is necessary to choose the basis functions {fi{z)} for the elliptic functions in order to further investigate 
the properties of the solution. A simple choice of the basis is 

fi{z) = 1, /2(z) = z, fsiz) == , f4(z) = , (3.44) 

Z Zif Z ~r Z-y 

where 



y(z) := ^(;,2_p2)(^2_q2). (3.45) 

Apparently, these are linearly independent on C since they have different pole structures. The coefficients 
{ci(^)} are determined by the conditions (|3.29p(|3.30p which can be written as a linear equation 

M-c = p, (3.46) 

where c and p are 4-component vectors which are given by Ci(^) and the residues, respectively, and M is 
a 4 X 4 matrix constructed in terms of the values /i(±^). The determinant of M is 



{e - zif 



detAf = 7;;'' ^^^^ . (3.47) 



For a generic ^, this is non- vanishing if and only if z,y is non-zero and finite. In this case, {ci(^)} are 
determined uniquely. The explicit expression of {ci(^)} is given in Appendix [Bl 
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The linear equation p.46p becomes singular, and therefore there is no solution, if Zy — 0,(X). This 
condition can be written as 

VT e A, (3.48) 

which was encountered as the condition for which the uniqueness of the solution may fail. 

Recall that the parameter ly was defined as p.20p . To realize a real n, v must take the following 
values, 

f (ill (0<n<2), , , 

V e { ^^' ^ ^ - '' 3.49 

[ 1 + iM (71 > 2). 

The sign of Im(i^) for n > 2 is irrelevant due to the uniqueness of the solution. 

The singularity at ;/ = 1 (n = 2) is anticipated since this is when the theory becomes degenerate. 
This singularity is special compared with the other ones since it is singular for any value of t. One might 
expect that a limit v —^ \ would be well-defined. In fact, a regular behavior in the 't Hooft coupling 
in this limit was observed in [18]. However, this is not the case at the level of the resolvents, as one 
can check using the explicit form of the coefScients shown in Appendix [BJ The individual coefficients are 
apparently divergent in the limit v ^^ 1. At the level of the 't Hooft couplings, there is a possibility of a 
cancellation of singularities among the coefficients. Since fi{z) are linearly independent, the resolvent is 
not well-defined in the limit. 

If n > 2, then there are many values of r for which the solution does not exist. Since the solution is 
well-defined for a generic value of r, the theory itself should make sense, and these singularities must be 
meaningful. In fact, these singularities appear as divergences of the 't Hooft couplings, as will be shown 
in subsection 13.41 Another viewpoint on these singularities will be given in section S) 

It has turned out that the solution of p.2ip is uniquely determined. In fact, this also implies that the 
resulting solution uj{z,^) has the following inversion property, 

oj{z-\o = -r'^{z,r'), (3.50) 

which, combined with 

r{z-\o = -r'nz,r'), (3.51) 

implies the inversion property (|2.23p . To show this, note that u{z) satisfies 

u{z-^) = u{z) + ^. (3.52) 

See Appendix O This implies G{z~^) = G{u{z) + |). Since G{u + |) also satisfies p.2ip . this can be 
written as a product of G(u) with an elliptic function R'{u). Due to the same argument, one finds that 
u!g{z^^,^) is given as 

4(^"\0 = ■^^^G{z)[c[iOfi{z) + c',iOf2{z) + c',iOf3{z) + c'^iOUiz)]- (3.53) 



Noticing a trivial identity 

1 111 



one finds that the residues of uj{z~-^,(,) at z ~ ±^~^ are related to the residues of U!{z, £,) as 



(3.54) 



Res±j-ia;(z-\^) = -C^Res±^uj{z,^). (3.55) 
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This implies that the coefficients {c^(C)} are related to {ci{^)} as 

cm = -e-'Q(r'), (3.56) 

which then implies p.SOp . 

3.3 Chern-Simons-adjoint theory revisited 

In [18j . Af = 3 Chern-Simons-adjoint theories were discussed. These are theories in C a member of which 
corresponds to a diagram Fad consisting of only a single node with an arbitrary number of ad edges. 
As explained in subsection 12.21 the solution in [TB] can be reproduced from the solution obtained in this 
section by imposing a suitable constraint. 

We have considered a theory F consisting of two nodes and a number of bf edges. The constraint 
requires that the eigenvalue distributions of two nodes of F are the same, making bf become equivalent 
to half the number of ad. In terms of the resolvent, the constraint is 

v\z) = v^{-z). (3.57) 

This implies that the vector-valued function w(z,^) must satisfy 



W(2,0 = ^i-Z,0 



This is realized by the pull-back of il{u,£,) satisfying 



1 

1 



(3.58) 



Using p.22p . this can be written as 



n^ {-u+^] = 02(u). (3.59) 



^s("+9) = -e^'^^Un)- (3.60) 



If ri5(M) satisfies this equation and the second equation in p.2ip . then the remaining equation is auto- 
matically satisfied. ^|(u) is then determined by the equation p.22p . 
The function G{u) satisfies 

g(u+Ij " e'^*'^G'(w). (3.61) 

If r2^(ii) is written as in p.24p . then one finds that R{u) is not just an elliptic function but it satisfies 

rL + ^J = -R{u). (3.62) 



As a function of z, this implies 



ciiO - C2(C)z + csiOhiz) + c,i0f3{z) 

- [ci(0 + C2{0z + csiOfziz) + cmUiz)] . (3.63) 
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Therefore, the solutfon of a Chern-Simons-adjoint theory is written by the coefficients {ci(^)} satisfying 

ciiO = 0, C3(e) = -C4(e). (3.64) 

The remaining coefficients are determined by requiring that ^^(2,^) has the appropriate pole structure. 
Since aj^(z,^) is now given by aj^(z,^), there are only two conditions. 

For the Chern-Simons-adjoint theory, one has to have ki = K2 since there is only one V{N) gauge 
group factor. The explicit formulas for the coefficients shown in Appendix [BJ show that this is compatible 
with p.64p . Note that the explicit functional form given in [T^ is different from the one obtained in this 
paper. This is simply because the choice of G(z) in this paper is different from the one in J18j . Since the 
solution is proved to be unique, these two expressions define the same function. 

3.4 't Hooft couplings 

Recall that the 't Hooft couplings are given in terms of resolvents as 

r = -v''(0). (3.65) 

This can be written in terms of cj(z,^) as in ()2.42|) . It is more convenient to consider the transformed 
coupling ts := {t^,t'^)S. For example, tg is given as 



t's = - d^ 



where 



ic 
1 



.. +*^— T^)777^ + f^5(0,e) 



^5(0,0 = -^G{0) 



cos ^ sin ^ J ^(l- 






(3.66) 



(3.67) 



e 

One can check that the integral is well-defined. Indeed, the integrand behaves as ^^^ in the limit ^ — > 00, 
while it is finite in the limit ^ ^- 0. 

The 't Hooft coupling ts depends on the positions p, q of the branch points. In fact, since pq — I 
is assumed, ts is a vector- valued function of p, say, and the domain of ts(j)) is the unit disk D. The 
behavior of ts{p) for a generic p would be quite complicated. We are interested in the possible divergent 
behavior of ts{p). In the following, we investigate not just the positions where is(p) diverges but also 
how it diverges. 

In ABJM theory and other Chern-Simons-matter theories with known gravity duals, interesting be- 
haviors are observed when the 't Hooft couplings diverge. Recall that imaginary values of {f} correspond 
to the physical values for the original Chern-Simons-matter theories. Therefore, a divergent behavior of 
tsip) would be interesting in the context of the original Chern-Simons-matter theories if both t^{p),t'^{p) 
have divergent imaginary parts while the real parts could be tuned to vanish in the limit. 

Apparently, ts{p) diverges when z^ = 0, 00. These are simple poles of ts{p) in D. It can be shown 
that the condition for Zi, = is 

lyr e A. (3.68) 
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Obviously, there is no solution for < n < 2 since for this range ^ < v < I. For the case n > 2, there 
may be many solutions. In fact, if r is purely imaginary, which is the case when p is real, then there are 
infinitely many solutions. This implies that ts{p) has infinitely many simple poles on the real axis. 

The fact that there are infinitely many simple poles for ts{p) implies that there are infinitely many 
values oi p for which ts{p) may have a prescribed large value. That is, for a given value of ts, there are 
infinitely many saddle-points for the original partition function p.2p . In this case, one has to choose one 
of the saddle-points having the smallest free energy. We will not discuss further about the criterion for 
choosing the dominant saddle-point. Naively, it could be expected that the minimum of the free energy is 
of order k^ (recall that Na ex fc in our planar limit) which would come from a saddle point corresponding 
to a value of p away from the origin. Similar results were found in |15| . 

The other possibilities for the divergence of ts{p) are also simple. Since ts{p) is given in terms of 
contour integrals, divergences appear when a singularity of the integrand approaches one of the endpoints 
of the contour, or when the contour is pinched by two singularities of the integrand. 

An endpoint singularity appears when the branch point at z = p in Ra{z) approaches the origin. Note 
that the other endpoint, at infinity, does not provide any divergence, and G{z) does not have singularities 
which can approach the origin. If p is very close to the origin, then a large contribution comes from the 
term 

K-i 4- K-^ 1 1 

(3.69) 



1 

1 



y(C) 



2 cos 2^ e 
in the integrand. The integral of this term gives a contribution 

which is divergent in the limit p — > 0. It turns out that t'g{p) also diverges in the same way. The original 
couplings {f°} then behave as 

The divergent part of the right-hand side is always real. To obtain large imaginary values for i°, one 
has to continue p along a path which winds around the origin many times, as was found in |18| for 
Chern-Simons-adjoint theories. The vanishing real part will be ensured only if \p\ = 0(1). 

Note that the limit p — >■ above is assumed not to be taken along the real axis for the cases n > 2, 
since there are poles. If one takes a limit such that 

T = Tl+iT2, |ti| = C < -, T2 ^> -l-OO (3.72) 

for a constant c, then the distance between vt and the lattice points in A is bounded from below, and 
therefore z"^ can be regarded as a smooth function oi p. 

The pinch singularity may appear when p approaches either ±i or — 1. In any cases, the resulting 
branch cuts, and therefore the eigenvalue distributions, have finite lengths. This is an uninteresting 
situation in the following sense. In ABJM theory, the expectation value of the BPS Wilson loop grows 
exponentially as the 't Hooft coupling grows. This is quite naturally explained from the dual string point 
of view [H] [10] [21] ■ On the other hand, if a Chern-Simons-matter theory has a finite size branch cut, 
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then the expectation vahie of the BPS Wilson loop corresponding to the cut must have a bounded value. 
The expectation value of a Wilson loop, say the one for U(iVa) gauge group factor, is given as 

where {zf } are the saddle-point values. The absolute value of this is bounded as 

\{W'^)\ < \q\ (3.74) 

since ztq are the eigenvalues with the largest absolute value. Therefore, if the branch cuts are of finite 
size even in the large 't Hooft coupling limit, then the magnitudes of the quantities (VF") are bounded in 
the limit. It seems to be difficult to explain this behavior by a smooth worldsheet in a classical geometry. 

In summary, the divergent behavior of {f^} as functions of p has a rather simple pattern. One might 
wonder whether the properties for the theories with Ug = 2 found so far would be generic ones for the 
theories in Cq. We will investigate a generic non-degenerate theory in Cq in section |4] and again find a 
similar behavior of the 't Hooft couplings and the Wilson loops. It might suggest that any non-degenerate 
theory in Co would not have simple gravity duals with the standard dictionary between the field theory 
and the gravity. 



3.5 Adding fundamental matters 

Now, let us consider the effect of adding f to the theory discussed so far. The homogeneous equations 
p.40p are the same, and therefore the functions {ujg{z,^)} can be written as 



4(z,e) = RHz))G{u{z)), u;l{z,0 = -R{-uiz))G(^-u{z) + ^ 



(3.75) 



The difference from the analysis in the previous subsections appears in the form of R{u). Since the 
matters f add poles at z = ±1 to vf{z) in (|2.41|) . the elliptic function R{u{z)) should have the following 
form 



^("(^)) = 77^^Tr^^2T^i(^(^)) (3-76) 



1 

(z2 - l)~z2 - ^2) 

with another elliptic function i?i (u) . It turns out that the divisor of Ri (u) is constrained as 



(Riiu)) >-(2,)-(-z,)-3(oo)-3(55). (3.77) 

The Riemann-Roch theorem implies that such functions form an eight-dimensional linear space. The 
residues of the extra poles at z = ±1 provide four more conditions. Therefore, the solution is again 
completely determined by the data given so far. 

The basis functions for Ri{z) are fi{z), • • • , fiiz) given in (|3.44|) and 

Mz) = z', k{z) = y{z), h{z) = z\ fs{z) = zy{z). (3.78) 
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The existence of the solution depends on the determinant of an 8 x 8 matrix M whose components are 
the values of the basis functions at the poles. The determinant is 

New singularities would appear when z^ = ±1. Again, this singularity does not appear when < 71 < 2, 
while it would appear for some particular values of r when n > 2. 

The addition of f introduces two poles at z = ±1 to u!s{z, £,)■ The functional form of a;5(z, 5) is different 
from the one in the absence of f. However, the behavior of t°-{p) would not be drastically modified by 
the addition of f since i" is determined by the values ujg{z, ^) at z = which is away from the new poles. 
Note that, although not directly related to our cases, in the case of the flavored ABJM theory [SI] the 
effects of f appear, for example, in a modified coefficient in the exponent of the expectation value (W) of 
the BPS Wilson loop. However, the qualitative properties of {W) as a function of the 't Hooft coupling 
are kept intact. 



4 Non-degenerate theory F in Co 



In this section, we investigate the planar solution of a general non-degenerate theory F in Cg. As was 
shown in subsection 13.51 the solution of a theory F including f can be obtained from the solution of 
another theory without f through a minor modification. In the following, the diagram F considered is 
assumed to have only bf. 

Let us recall the technique used to solve a non-degenerate theory F with two nodes. The resolvents 
{v"'{z)} were obtained from the pull-back of functions {^"'{u)} defined on the u-plane. The relation 
between the coordinates z and u was given by the Abel-Jacobi map 



uiz) = / A. (4.1) 

Jp 

This actually defines a map from the elliptic curve E to its Jacobian J{E) which is again a torus. The 
equations p.l7p suggest that i}{z) would be given by a section of a vector bundle on J{E). The vector 
bundle is specified by the monodromy matrices / and M2M1 for the A-cycle and the _B-cycle, respectively. 
By diagonalizing the monodromy matrices, the problem is reduced to the one for a pair of line bundles 
on J{E). Since the theta functions are sections of some line bundles on J{E) [JD], it is natural that VL{z) 
can be written in terms of the theta functions. 

Now, consider a theory F with many nodes. The equations (|2.40p . supplemented by the conditions 
for the absence of branch cuts, can be written as 

Uj{yt,S,) ^ ^{Va^OMa, ya&[Pa,qal (4.2) 
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where we have assumed PaQa — 1, and 



M„ 



1 n"'"-^ 
-1 

^a,a+l I 



(4.3) 



Note that each Ma is diagonaUzablc and satisfy M^ = /. 

One may consider, instead of the eUiptic curve E, a higher-genus Riemann surface S defined by 



y 



Y[(x-pa){x-qa)- 



(4.4) 



The genus of S is 7ig — 1. The Abel-Jacobi map can be defined for higher- genus Riemann surfaces [40] . 
Let A' (i = 1, ■ 



, rig — 1) be a set of independent holomorphic 1-forms on E. Then, a set of integrals 



u\z) 



X' 



(4.5) 



where p is a point on C, defines a map from S to its Jacobian J(S) which is a torus of dimension Ug — 1. 
This map is an embedding of E into J(E). The A-cycles and the B-cycles of E are mapped to non-trivial 
cycles of J(E). 

A possible generalization of the strategy used for the case Ug = 2 would be to construct a vector 
bundle on J(E) and then to take the pull-back of a section of the vector bundle. The vector bundle 
should be characterized by some monodromics around the non-trivial cycles on J(E). In analogy with 
the case Ug = 2, the monodromics for the cycles of J(E) corresponding to the A-cycles of E would be 
trivial, while those corresponding to the i3-cycles would be given by matrices like MaMa-i- In this way, 
the equations (j4.2p would define a homomorphism 



p : ^i(J(E)) ^ GL{ng,C), 



(4.6) 



that is, a representation of the fundamental group tti ( J(E)). The problem is that the matrices {MaMa-i} 
do not commute among them. Since 7ri( J(E)) is Abelian, such a representation p cannot exist. Therefore, 
a generalization along this line is not successful. 

It would be still possible to regard aj(z,^) as a section of a vector bundle on E, not on J(E), with 
monodromics obtained from Ma- The fundamental group tti (E) is a non- Abelian group which is generated 
by 2(ng — 1) letters ai, /3i, • • • , an^-i, Pug-i subject to a single constraint 



n a.l3^a-'|3^' = 1- 



(4.7) 



Note that from this description it is easy to notice that the fundamental group of a torus is Abelian. 
As mentioned above, there must exist a homomorphism from 7ri(E) to GL(ng,C). A map defined by 
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assigning the unit matrix to each ai and a suitable product of Ala to each /3i apparently preserves the 
relations (|4.7p . and therefore defines a homomorphisni. 

The desired vector bundle on S can be constructed as follows [41]. Let S be the universal cover 
of S. E can be regarded as a principal bundle P over E whose structure group is 7ri(I]) [421 . If an 
rig-dimensional representation p of tti(T,) on V^ = C"" is given, then one can construct the associated 
vector bundle P XpV which has the desired monodromy. The construction of the vector bundle is rather 
simple, however, the description of sections of the vector bundle does not seem to be tractable. 



4.1 Riemann-Hilbert problem 

It turns out that there is a more convenient reformulation. The equations (|4.2p can be regarded as an 
implication of the existence of monodromies for a;(z,^) around the points {PaiQa}- The monodromy 
matrices at both z = Pa,qa are Ma since M~^ = Ala holds. Now the problem of solving the equations 
(|4.2p becomes constructing a vector bundle on a punctured sphere 

V := Pi\{pi,gi,...,p„^,<7„j, (4.8) 

whose (multi- valued) section a;(z) has the monodromies around the punctures prescribed above. 

As in the case of the higher-genus Riemann surface discussed above, the construction of the vector 
bundle is rather straightforward. The advantage of this reformulation is the relation to a well-studied 
mathematical problem, known as the Riemann-Hilbert problem, in its original form. This relation then 
provides us with a rather detailed description of sections of the vector bundle. A review of the Riemann- 
Hilbert problem can be found, for example, in (43] and |44], the latter of which also includes recent 
developments. A pedagogical introduction to the basics on analytic differential equations can be found 



The Riemann-Hilbert problem in the original form discusses a relation between differential equations 

and monodromies exhibited by the solutions of them. Let S^ be the set of n-dimensional Fuchsian 

systems, a system of n first-order differential equations with Fuchsian singularities only. Let pi {i = 

1, • • • , to) be points on C Then, a Fuchsian system in iS„ whose Fuchsian singularities are at {pi} is of 

the form 

ay >r- A, 

-T- = V y (4-9) 

dz "^ z~ pi 

?— 1 

where y{z) is an n-vector-valued function and {Ai} are constant n x n matrices. The Fuchsian system is 
non-singular at infinity if and only if 

rn 

5]A, = (4.10) 

1=1 
is satisfied. It is always assumed below that the point at infinity is a non-singular point. 

There exist n independent solutions {y\z)} (Z = 1, • • • , ri) of (|4.9p which are holomorphic on V. As is 

familiar for the hypergeometric function, these solutions exhibit a monodromy when they are analytically 

continued around a singularity. It is convenient to consider the fundamental matrix 

Y{z) := (yi(z),...,2/"(z)) (4.11) 
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where the sokitions {y'(z)} arc regarded as column vectors. Note that Y{z) is non-degenerate over D. 
When Y(z) is continued counterclockwise around a singularity z = p^, it transforms as 

Y{z) -> Yiz)M, (4.12) 

with a non-degenerate constant matrix Mi. Typically, logAf^ is proportional to Ai up to a similarity 
transformation. 

Let 7i denote a closed path, starting and ending at a chosen point, and encircling only one singularity 
z = Pi once and counterclockwise. Suppose that the singularities are numbered such that the composition 
of the paths 

71 072 o • •• o7,„ (4.13) 

is homotopically trivial on V. Then, the monodromy matrices must satisfy 

M^---M2Mi - /. (4.14) 

This construction defines a map 

0RH : 5f -> 5f , (4.15) 

where S^^ consists of sets of pairs {pi,Mi) satisfying the condition (|4.14[) . An element {{pi,Mi)} of 
S^^ is called a monodromy data. When a monodromy data {{pi,Mi)} has a non-empty inverse image 
(j)j^-^{{{pi, Mi)}), it is said that the Riemann-Hilbert problem for this monodromy data has the positive 
solution. 

The Riemann-Hilbert problem was originally concerned with the surjectivity of (/)rh- At first, it was 
believed that 0rh is always surjective. However, various counterexamples have been found |43j . Up to 
now, various sufficient conditions for the positive solvability of the Riemann-Hilbert problem have been 
found. One simple sufficient condition is the following |46j : if at least one of the monodromy matrix in a 
given monodromy data is diagonalizable, then the Riemann-Hilbert problem for this case has the positive 
solution. Another sufficient condition for the positive solvability of the Riemann-Hilbert problem is the 
irreducibility of the monodromy data [17] [33 • 

These theorems can be proved by relating the problem of the existence of a Fuchsian system to the 
problem of the existence of a holomorphic vector bundle on P^ with a meromorphic connection. It can 
be shown |45j that, for any monodromy data, one can construct a holomorphic vector bundle on P^ with 
a meromorphic connection V such that a horizontal section s, a section satisfying Vs = 0, exhibits the 
prescribed monodromy. The meromorphic connection can be also chosen such that it has only simple 
poles. If the holomorphic vector bundle is trivial, the equation Vs = in terms of a coordinate is 
equivalent to a Fuchsian system (|4.9p . 



When the Riemann-Hilbert problem has the positive solution for a given monodromy data, one may 
be able to obtain a matrix- valued function exhibiting the prescribed monodromy. This function is nothing 
but the fundamental matrix Yq{z) of the corresponding Fuchsian system. 

The structure of Yq{z) can be described explicitly [33]. Let Ei be a matrix related to the monodromy 
matrix Mi as 

M, = S-^eyi^{2T:iE,)S., (4.16) 
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with a suitable non-degenerate matrix Si. The matrix Ei is chosen such that it is of the Jordan normal 
form, and its eigenvalues {/3'} satisfy 

< Re(/3,^) < 1. (4.17) 

Define a matrix-valued function z^^ for a matrix M as 

z'^' := exp(A/logz). (4.18) 

In terms of the matrices Ei, Si, the local form of Yo{z) around z = pi can be written as 

Yo{z) = Uiz){z^p,)^'iz~p,f'S, (4.19) 

where U{z) is a matrix- valued function which is holomorphic and non-degenerate at z = pi, and A^ is a 
diagonal matrix whose components are integers. The matrix A; has to be chosen appropriately at each 
singularity such that Yo{z) is a section of a trivial holomorphic vector bundle. 

It can be shown [49] that there exists a matrix- valued function r(z) which is meromorphic on P'^ and 
holomorphic and non-degenerate on 'D\{oo} such that the product 

Y{z) := r{z)Yo{z) (4.20) 

has the local form (|4.19p with A^ absent at all singularities. Moreover, the function Y{z) has the following 
simple form near infinity 

Y{z) = diag(z=S--- ,z'=")T/(z) (4.21) 

where {ci} are integers and V{z) is a matrix-valued function which is holomorphic and non-degenerate 
at infinity. 

The integers {ci} have the following geometric meaning [33]. Let "D be the universal cover of V. 
As mentioned before, this can be regarded as a principal bundle on V whose structure group is the 
fundamental group 7ri(I'). One may construct a holomorphic GL(ri,, C) bundle P associated to the 
principal bundle V. The representation p : 7ri(2?) — > GL(n,C) necessary to construct P is defined in 
terms of the monodromy matrices {Mi}. A section of this bundle defines a holomorphic connection V on 
P. It is possible to extend P to a holomorphic GL{n, C) bundle on P^ with a meromorphic connection 
extending V. Each extension is specified by a set of matrices {Si,Ai}. The non-degenerate matrices 
{Si} are chosen such that the similarity transformation by Si makes Mi upper-triangular. The other 
matrices {Ai} are diagonal matrices with integer components. The canonical extension P° is defined as 
an extension of P for which all A^ are zero. It is known that P° does not depend on {Si}. Let F° be 
the holomorphic vector bundle associated to P° . The Birkhoff-Grothendieck theorem implies that any 
holomorphic vector bundle on P^ is holomorphically equivalent to a direct sum of line bundles. For the 
holomorphic vector bundle F° , one finds 

F° = C'(-ci)e---®C'(-c„). (4.22) 

The integers {c;} are known as the splitting type of F° . 
The sum of c; is determined by Ei as [50] 

n rn 

J2ci ^ Y^tiiE,). (4.23) 

i=l 4=1 
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Therefore, this sum can be determined solely from the monodromy data. An algorithm of calculating 
each ci is available |49j once one knows the Fuchsian system explicitly from which one started. 



4.2 The solution 

Let us return to the matrix model problem. The equations (|4.2[) imply that the monodromy data relevant 
to our case is A4 := {{pa, Ma), {qa, Ma)}. This monodromy data is consistent since 

l[M{pa)M{qa) = I (4.24) 

i=l 

is trivially satisfied. As mentioned below (|4.3p . any Ala is diagonalizable. Therefore, the corresponding 
Riemann-Hilbert problem has the positive solution, implying the existence of a matrix-valued function 
Y{z) with the properties described in the previous subsection. It would be interesting to notice that the 
condition (j2.43p is equivalent to the one for the reducibility of the monodromy. See Appendix [C] 
Each matrix Ea can be chosen to be the following diagonal matrix 

Ea = diagfio,--- ,oV (4.25) 



^2 
This implies that the function {z ~ pa)^"^ is also diagonal, 

(z - Pa)""-" = diag {V^^Va, 1, • • • , 1) ■ (4.26) 

The formula (|4.23p implies 

Ug 

Y.CI = ng. (4.27) 

1=1 

To determine the values of {q}, let us first discuss the uniqueness of the solution of ()2.33|) . Suppose 
that there are two solutions v{z) and v'{z) written in the vector notation. The difference 5v{z) := 
v{z) — v'{z) satisfies the homogeneous equations (|4.2[) . Therefore, the function 5v{z) has monodromies 
described by the monodromy data Ai. Then a product 

r{z) := 5v{z)Y{z)~^ (4.28) 

has no monodromies, implying that the components of r{z) are meromorphic functions on P^. Since 
both v{z) and v'{z) are assumed to be finite at the points {pa,Qa\, so is 5v{z). The determinant of Y{z) 
behaves like ^Jz — pa or ^J z — qa near the singularities, and therefore, Y{z)~^ cannot become singular 
enough to create poles in r{z). It is concluded that r(z) must be holomorphic on C 
The behavior of r(z) near infinity is 

r{z) = {Sv\z)z-''\--- ,Sv'^<'{z)z-'="'j), dv{z) := Sv{z)Viz)~\ (4.29) 

If all {c/} are positive, then this implies that r{z) vanishes at infinity. Since r{z) turns out to be 
holomorphic on P"'^ and vanishes at infinity, r(z) must vanish, implying the uniqueness of the solution of 
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(j2.33p . On the other hand, if some of {q} are zero or negative, then there are ambiguities of the solution 
which cannot be fixed by the conditions given so far. 

We found in section [3] that the solution is unique when ng = 2. This implies that in this case ci, C2 > 0. 
Since ci + C2 = 2 must be satisfied, one finds ci = C2 = 1. For all ng > 2, there exists a theory in Co for 
which the existence and the uniqueness of the solution can be shown [51j . Therefore, for those theories, 
one must have 

ci = ••• - c„^ = 1. (4.30) 

Since {q} are integers, they would not change under continuous deformations of the holomorphic vector 
bundle F°. If {n°-^} are regarded as complex numbers, then one may be able to consider continuous 
deformations of {Ma} which induce corresponding deformations of F°. Note that the Riemann-Hilbert 
problem is still valid for such complex values of {n"''}. Since any possible {Ma} are connected to the 
monodromy matrices for a theory discussed in [51], the splitting type {ci} would be (|4.30p in general. 
This implies the uniqueness of the solution. 

Note that this conclusion would be invalid if there would exist a domain wall in the complexified 
parameter space {n"-''} across which the values of {ci} jump. It would be interesting to clarify this issue. 

If the values of {q} are given as (j4.30p . then it turns out that Y{z) has a nice inversion property. 
To show this, consider Y{z~^). This has the same monodromy property as Y{z), due to the relation 
PaQa — 1- This implies that a product Y{z~^)Y{z)~^ has no monodromy. This product has a simple pole 
at the origin and a simple zero at infinity. This implies 

y(z-i)y(z)-i = z-^c (4.31) 

with a constant matrix C. Since the equality should hold at z = 1, it is concluded that 

Y{z-^) = z-^Y{z) (4.32) 

holds. At first sight, this would not be valid at z = —1. This is actually not that case since Y{z) has 
branch cuts. As one takes z = e*^ and changes the value of 6 from to tt, one finds that z and z^^ pass 
through those cuts, resulting in the extra minus sign. 

The function Y{z) does not change the monodromy if it is multiplied by a constant matrix from the 
left. Using this degree of freedom, we normalize Y{z) such that Y{0) is the unit matrix. 

Now, let us solve the equations (|4.2|) . Suppose that there is a solution a;(z,^). Then, the product 
of this solution with Y{z)~^ multiplied from the right has no monodromy. Therefore, they must be a 
vector- valued function whose components are meromorphic functions on P^. In fact, they are rational 
functions. Let r{z,^) denote this vector-valued function. Then a;(z,^) can be written as 

uj{z,0 = riz,OYiz). (4.33) 

The function r(z,^) is determined by requiring that lo{z,^) has the appropriate pole structure, as in 
section[3l Recall that the function Y{z) is defined such that Y{z) is finite at the singularities z = pa, Qa- 
Since we have assumed that io{z,^) is finite at z = pa,qa, it turns out that ^(z,^) has poles only at 
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z — ±^. In addition, r(z,^) should behave as z ^ near infinity since Y{z) diverges as z while 0-1(2;,^) was 
assumed to be finite. Therefore, r{z^i^ must have the following form 

'^(^'O = ^^^[«(e)~^ + &(C)] (4.34) 

where a(^), 6(^) are z-independent row vectors. Let the required residues of a;(z, ^) at z = ±^ be denoted 
as 

p± := Res±jw(z,e), (4.35) 

The vectors a(^), 6(^), c(^) are determined by the equations 

a(e)e + fc(0 = 2p^iY{i)-\ (4.36) 

<i)^-Ki) = 2p^iY{-i)-\ (4.37) 

These equations can be solved easily, and the solution is 

a(0 = p^Y{i)-^+P-Y{-i)-\ (4.38) 

Ki) = i{p+Y(£)-^-p^Y{-i)-^\. (4.39) 



Therefore, the solution can be written as 



w(z,0 



^ p+y(e)-i + ^p_y(-0-^ 



As was argued above, this is the unique solution of the equations (|4.2 



y(z). (4.40) 



4.3 Local theory of Fuchsian systems and 't Hooft couplings 

It will be quite difficult, unless impossible, to write down the explicit form of aj(z, ^) in terms of well-known 
functions. However, the solution (|4.40p obtained above is not useless. It was found in subsection l3.4l that. 
as long as we are interested in the divergent behavior of the 't Hooft couplings, the only information 
necessary to extract it is the local behavior of the solution around the branch points. Fortunately, the 
local structure of ^(z), which is the only non-trivial function in the solution (|4.40p . can be determined 
as (|4.19p with Aa = based on the local theory of Fuchsian systems [33] . 

The 't Hooft coupling t = — 1'(0) is given as 



t = - \ de[r(0, 0+^(0,0], (4.41) 

Jc 

where 

^(0, ^) = \ V-p^^iiY^ + p-y{-ir^\ ■ (4.42) 

For the case rig = 2, the divergent behavior of interest appeared when one of the branch point approaches 
the origin, one of the endpoint of the integration contour C. 
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Let us investigate the behavior of t when a branch point, say Pa, approaches the origin. The deter- 
minant of Y{^) near the point S, — Pa ^s 



dcty(e) = dctUiOVl^Va (4.43) 

since in our case tr(£'a) — i. Therefore, if Pa is very close to the origin, then the behavior of y{£,)~^ foi' 
smah ^ should be 

Y{0-' - \[t^^^ (4.44) 

due to the normalization Y{Q) = I. This implies that the divergent behavior for a non-degenerate theory 
is quite similar to the one observed for the theory with Ug = 2. In fact, the 't Hooft coupling t diverges 
as 

t {p+-p^)\ogpa (4.45) 

as Pa approaches the origin. 

In general, it would be possible that several branch points approach to the origin simultaneously. Even 
in such cases, the divergent behavior is similar. One may find that the integral (j4.4ip can be estimated 
roughly as 

t - -j^^{-p+ + p-) (4.46) 

where e is determined by the parameters {pa}- Therefore, in any case the divergence is logarithmic, and 
in particular, only the real part of t may diverge. 



4.4 Isomonodromic deformations 

It has been shown that the solution for a non-degenerate theory with Ug > 2 can be obtained, and the 
solution exhibits a similar divergent behavior of the 't Hooft coupling t as observed in section |3] for the 
theories with Ug = 2. There was another observation for the theories with Ug ~ 2. That is, the solution 
may not exist for a particular choice of the parameters. This might seem to be contradicting with the 
conclusion obtained in this section. Indeed, we found the solution of the equations (|4.2p for any n"-^ and 
for any positions of the branch points. The existence of the solution was based on the positive solvability 
of the corresponding Riemann-Hilbert problem which was rigorously proved |46j . 

The resolution comes from noticing that the solution in section |3] depends analytically on the positions 
of the branch points. The theorem in |46| claims that, for any configuration of the singular points, there 
exists a Fuchsian system reproducing the prescribed monodromy. However, this does not claim that the 
family of the Fuchsian systems is analytically parametrized by the positions of the singularities. In fact, it 
is known |52j j 53 j [54] that, although some of the Fuchsian systems form a continuous family parametrized 
by the positions, the parameter space is not the entire C™, where m is the number of singularities of 
the Fuchsian system. For a particular choice the positions there is no Fuchsian system in the continuous 
family which reproduces the prescribed monodromy. The divergence of the 't Hooft coupling t should be 
related to this phenomenon. 
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The relevant notion here is the isomonodromic deformations 
the foUowing family of Fuchsian systems 



dy 

dz 



Ma) 



Z — J y — 



-y, 



of Fuchsian systems. One considers 



(4.47) 



where the matrices {Ai{a)} depend on the positions {a;} of the singularities. The dependence of A,;(a) 
on {fli} is chosen such that the monodromy realized by the solution of (|4.47p does not depend on {ui}. A 
suitable change of Ai{a) according to the continuous change of a^ preserving the monodromy is called an 
isomonodromic deformation. The family of the Fuchsian systems (|4.47p should be the appropriate one 
for discussing the solution of the equations (|4.2p since we would like to discuss physical quantities, for 
example the 't Hooft coupling t, as a function of the positions of the branch points. 

The matrices {Aj(a)} for an isomonodromic deformation are known to satisfy the Schlesinger equations 



dA,{a) 



^[Aia),A,ia)],^^^ 



a,: 



a 



aj). 



(4.48) 
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It was shown [52] [53] [54] that the solution of this equation has movable singularities. A movable singularity 
is a singularity whose position depends on the initial condition for the solution, and therefore, the existence 
of such singularities is not easily anticipated from the equation. If a given set {ai} of values corresponds 
to the position of a movable singularity, then this implies that the Fuchsian system exhibiting a prescribed 
monodromy at the singularities {a^} does not exist within the continuous family of Fuchsian systems. 

For the case Ug = 2, the above discussion on movable singularities can be made more explicit. It is 
known (see e.g. [5^) that in this case the Schlesinger equations are reduced to the Painleve VI equation 
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(4.49) 



Note that the relation between matrices in the Fuchsian system and the parameters a, /3,7,(5 in (|4.49p 
was given, for example in |57j . for the case when the monodromy matrices belong to SL(2,C), but in 
our case this is not since detAf^ = —1. Actually, this is not a big problem since any monodromy in 
GL(n, C) can be decomposed into a one-dimensional monodromy and a monodromy in SL(n, C), and the 
one-dimensional monodromy problem can be easily solved. In the case when all monodromy matrices 
belong to SL(2, C), the parameters a, (3, 7, 5 are determined by the eigenvalues iia^ of the matrices {Ai} 
in the Fuchsian system as 



a 



(ai - 1)^ 



/? 



7 



a% 



5 



(4.50) 



2 ' ' 2 ' 2 2 2 

The matrix Ai is equal to a matrix E'i, which is related to the monodromy matrix as in (|4.16p . up to a 
similarity transformation. In our case, one finds ai = j, implying 



(3 



1 



1 



1 



S=l 



(4.51) 



Quite remarkably, even though the Painleve VI equation ()4.49p is highly non-linear, a two-parameter 
family of solutions of this equation is known |58| for this choice of parameters. It is written in terms of 
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theta functions. The explicit form is given in the foUowing parametric form 



2/(t) 



t?i"(0) 1 



37r2^4(0)X(0) 3 



^ , ^3(0)^ 



^4(0)^ 



i?'/'(i/)i9i(j/) - 2&[{v)d\{v) ^ ^■Kicx[d'{{v)^x{v) - d\[vf) 



(4.53) 



where the modulus t of the thcta function is the modulus of the elliptic curve defined by 

4 

y" = n^^-"^)' (4.54) 

and V here is related to the parameters C\ , C2 as 

V — c\T + C2- (4.55) 

It is apparent that y{x) has a pole at x for which 

C1T + C2 e Z + Zr. (4.56) 

This is quite similar to the condition p.48p for which there is no solution of the equations p.40p . and 
for which the 't Hooft couplings diverge. Exactly speaking, the modulus r here corresponds to r^^ in 
section [3J but this does not make any essential modification. 



5 Discussion 



We have investigated the planar solution of various Chern-Simons-mattcr theories. Under the assumption 
of the symmetric distributions of the eigenvalues, one can reduce the problem of solving a generic Chern- 
Simons-matter theory to the one of solving the corresponding theory including only fundamental and 
bi-fundamental hypermultiplets. Due to this planar relations, it was shown that many of the theories in 
the family C defined in section[2]may admit a systematic analysis of the planar resolvents. As a result, we 
obtained an integral representation of the resolvents whose integrands are written in terms of the solution 
of a Fuchsian system, a set of first-order differential equations, as long as the theory is non-degenerate, 
that is, the condition (|2.43p is not satisfied. Local properties of the solution of Fuchsian systems have 
been studied so far. The results can be used to clarify the properties of physical observables of the 
Chern-Simons-matter theories. It was found that the properties of the non-degenerate theories are quite 
different from those of ABJM theory For example, the vevs of Wilson loops do not grow exponentially 
as the (physical) 't Hooft couplings grow. These theories are rather similar to the Chern-Simons-adjoint 
theories investigated in [18]. For the case when the gauge group is the product of two U(7Vq) factors, the 
integrands for the resolvents can be written more explicitly in terms of theta functions, and the properties 
of the solution can be analyzed in detail. 
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Figure 7: The Chern- Simons-matter theory corresponding to this diagram can be solved by the technique 
in this paper. This is reahzed by a hncar array of D5-branes and NS5-brancs with D3-brancs suspended 
among them. 

An interesting resuh obtained in this paper is the condition (|2.43p . According to the resuhs in section 
m it would be expected that an interesting large 't Hooft coupling limit exists only for degenerate theories. 
Note that the Chern-Simons-matter theories satisfying the "no long-range force" condition discussed in 
[30] also satisfy (|2.43p . From the point of view of the matrix model, (|2.43p is a rather technical condition 
which shows the inability to apply a particular technique developed in this paper. It would be very 
interesting to clarify the meaning of this condition in the context of quantum field theory. Since this 
condition seems to be a necessary condition for Chcrn-Simons-mattcr theories to have simple gravity 
duals, the understanding of this condition would provide us with an insight on the mechanism of the 
emergence of a dual geometry in the context of AdS/CFT correspondence. The condition (I2.43|) should 
be also important to understand the large N behavior of the free energy which scales as A^ 2 for ABJM 
theory [llj in terms of quantum field theory. Since such a behavior is not universal for general Chern- 
Simons-matter theories, the explanation should take into account some effects of a particular choice of 
the matter contents. 

It would be also interesting to classify the solutions of ()2.43|) and to investigate the properties of 
the corresponding theories. It was shown in subsection 12.41 that the solutions are not exhausted by n"-'' 
corresponding to ADE afhne Dynkin diagrams. To analyze those degenerate theories, one has to invent a 
new technique for solving the saddle-point equations for those theories. Probably, some of those theories 
can be solved more completely by using the Fermi gas approach [T^. It would be interesting to specify 
the subset of those degenerate theories to which the Fermi gas approach is applicable. 

It is usually quite helpful to find a brane realization of a Chern-Simons-matter theory for understanding 
its properties. Of course, one should not expect that all the theories in the family C might have such 
a brane realization. Since many of the theories with known gravity duals satisfy (|2.43|) . one might 
anticipate that brane realization might not be available for non-degenerate theories. However, it would 
be interesting to notice that there exist some non-degenerate Chern-Simons-matter theories with brane 
realizations. An example is the theory corresponding to a linear diagram depicted in Figure [T] With the 
help of the knowledge on the brane dynamics, one may be able to investigate such a theory in more detail 

An interesting physical issue for Chern-Simons-matter theories is supersymmetry breaking. For simple 
theories, it is possible to investigate the condition for when the supersymmetry is spontaneously broken 
[59] [60] [61] [62] . Recently, it is argued that the free energy (or the partition function) can be regarded 
as a kind of an order parameter of supersymmetry breaking [63| . For example, for pure Chcrn-Simons 
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theory, the partition function [5] 

(sin—) (5.1) 

m— 1 

vanishes when the supersymmetry is broken, that is, N > k. It would be interesting to evaluate the free 
energy using the planar resolvent obtained in this paper, and systematically investigate singular behaviors 
of the free energy to check whether it is really related to supersymmetry breaking. 

The possible dual to a Chern- Simons-matter theory is not only a gravity theory, but it may be a 
higher-spin theory [M]. It would be interesting to find a criterion by which one can judge whether a 
given Chern-Simons-matter theory may be dual to a higher-spin theory. 
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A The Abel-Jacobi map 



Let Ao be the following holomorphic 1-form 

Ao := '^^ =. (A.l) 

on the elliptic curve E defined by 

y' = {x'-p'){x'-q'). (A.2) 

The integrals of Ao 

TTA := / Ao, ttb := / Aq (A.3) 

J A J B 

over the A-cycles and the i?-cycles in Figure [5] satisfy the second Riemann relation 

Re(7rA)Im(7rs) - Re(7rs)Im(7rA) > 0. (A.4) 

Usually, a normalized 1-form is used for which tta — 1 such that Im(7rs) > 0. In section [3l we instead 
employ another normalization 



A := — Ao, (A.5) 



for which Im(7r^) < 0. 
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The map 



u{z) := / A 
Jp 



defined in subsection 13.11 maps several points in the z-plane as follows, 



1 T 1 1 r 1 r 

2 - 2' 2'4' '~2'4 ^ 2 



(A.6) 



(A.7) 



and a half of the z-plane is mapped onto a parallelogram formed by the above points. The other half is 
mapped to another parallelogram. 

It is obvious from the definition that 



u(z) = —u{z) 
holds, where z is a point on the second Riemann sheet. It is also easy to show that 

P"^ 1 

u{-z) = A = -u{z) + -. 

Jp 
Combining these two relations, one obtains 

u{T) = u{-z)- -. 

The u-coordinate of z~^ should be given with care. Naively, 



A = f X+ f \ 

J p J q 



suggests 



u{z ^) = u{z) + 



(A.; 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 



A precise definition would require the choice of a contour, and also the determination of which Riemann 
sheet z~^ lies on. We instead define z~^ such that the above relation holds. Note that this implies 
1-1 = 1 etc. 



B Coefficients Ci{^) 



The coefficients {ci{^)} in subsection 13.11 are determined by the following linear equation 



where 



1 C MO hiO 

1 -e f3i~o /4(-e) 

1 C hi-O fsi-O 

1 -e /4(e) fsio 



CliO 
C2{0 

C3(e) 

C4(0 



2i^ 



H2G+ 

K2G_ 
— KlG-l- 



G± := e±^"''G(±0^^ 



(B.l) 



(B.2) 
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The solution is 



ci(0 

C2(0 

C3(e) 

C4(0 



4z^ sinTTZ/ y(^) 
«(«;i+K2) 1 



i 1 



{A-G+-A+G-), 



i 1 



i>i2C+ + ^iG-)iG+-G-), 
(kiC+ + k2C_)(G+-G_), 



where 



4z, 


aimivyiC) 


V'-i^i 


- 1 :i.^^ — 


A± 


= 2z^[z,y 


y{z,y) 


±ey(e)], 


B± 


= 2Ky(z 


.)±^ 


^.2/(0], 


c± 


= (C±^- 


.)(^^ 





(B.3) 
(B.4) 

(B.5) 
(B.6) 



(B.7) 
(B.8) 
(B.9) 



These coefhcients arc finite at ^ = 0. Rccalhng that G{z) has a simple pole at infinity, Ci{£,) are also 
finite at ^ = 00. This implies that the integral (j2.35p are well-defined. 



C Irreducibility of the monodromy 

The monodromy representation discussed in section |4] is the one defined by the matrices {Ma} given as 
(|4.3p . This is irreducible if there is no proper subspace of C"« which is fixed by all Mq. 
Assume that there exists a proper subspace V C C"« such that 



V-Ma C V. 



(C.l) 



Let m < Ug he the dimension of V. One may choose a basis of V. The components of the basis vectors 



form an m X Ug matrix 



^m ' ^m 5 



(C.2) 



The basis can be chosen such that all m x m minors are non-degenerate. In fact, one choice of such 
a basis corresponds to a point of the Grassmannian Gr^^n (C) contained in the intersection of all the 
standard coordinate patches, which is obviously non-empty. 

Let w be a vector in V. It is parametrized by m variables {t*} as 



(elf 



,c:^n 



(C.3) 



Let us fix a value of a. By assumption, vMa is in V . Therefore, there is another set of parameters {?'} 
such that 

vMa = (cif,---,c-''f). (C.4) 
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This and the exphcit form (|4.3p of Ma impUcs that, for & 7^ a, 

c\{f - n = 0. (C.5) 

One can choose m of these equations for which c\ is non-degenerate. These equation imphes t* = i' for 
all i. 

The remaining equations are 

^ n"\¥ - 2c,^f = 0. (C.6) 

b 

There exists a non-trivial solution for {t*} if and only if 

det [2(5'^'' - n''''] = 0. (C.7) 

Therefore, for a non-dcgencrate theory for which this is not satisfied, the monodromy representation is 
irreducible. 
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